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The origin and nature of ultrahigh-energy cosmic ray events, above the 
Greisen-Zatsepin-Kuzmin (GZK) cutoff energy, constitute a long-standing, 
unsolved mistery. Neutrinos are proposed candidates but their standard 
interactions with matter are too weak. In the context of a TeV-scale string 
theory, motivated by possible extra space dimensions, the neutrino- nucleon 
scattering is examined. Resonant string contributions increase substantially 
the standard model neutrino-nucleon cross section. Although they seem 
insufficient to explain the trans-GZK cosmic ray events, their effects might 
be detected in next experiments. 

1. The mistery of ultrahigh-energy cosmic rays 

Cosmic rays (CRs) were discovered in balloon experiments by V. Hess 
in 1911. They have always been a fruitful particle physics laboratory: 
positrons, muons, pions and many other particles were first observed there. 
CRs have a broad spectrum, ranging from 1 GeV to 10^*^ eV, including en- 
ergies far beyond the ones available at man-made accelerators. The flux of 
CRs decreases rapidly with the energy: e.g. 1 particle/m^/s at 10^^ eV and 
only 1 particle/km^ /year at about 10^^ eV. 

At high energies, a CR may hit an atmospheric nucleon producing a 
cascade of secondary particles, extensive air showers (EAS), that cover at 
ground an area with a diameter of tens of meters at 10^^ eV to several kilo- 
meters at 10'^'' eV. Incidentally, it is easier to observe the EAS at ground 
level for energies where direct detection of CRs at satellite or balloon ex- 
periments becomes impractical. The primary particle (mass and energy) 
can only be inferred from the parameters of the shower, such as number of 
electrons, muons or hadrons, shapes of lateral distributions, etc. 
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We will focus on the part of the spectrum above the Greisen-Zatsepin- 
Kuzmin cutoff energy [1] Eqzk ~ 5 x 10-^^ eV (see below). Unexpectedly, 
over seventy trans-GZK ultrahigh- energy cosmic ray events have been ob- 
served over the last forty years by five different experiments [2]. These 
events exhibit large-scale isotropy and small-scale anisotropy (a very un- 
likely pairing within the angular resolution occurs: three doublets and one 
triplet!). They constitute a puzzle in modern astrophysics. The main open 
questions are: 

• What is the primary particle? 

Protons with energies above the reaction threshold Eqzk lose part of 
their energy by scattering with cosmic microwave background (CMB) 
photons: 

P + 72.7K —>■ —>■ p + tt'^ {n + 7r+). 

Above a distance of Dqzk ~ 50 Mpc the energy of the proton gets 
below lO^'^ eV, regardless of its initial value, after a series of such col- 
lisions. This implies that the source cannot be too distant for trans- 
GZK events if they are originated by protons. Similar cutoffs exist 
at lower energies for other primaries: nuclei undergo photodisintegra- 
tion in CMB and IR radiations; electrons lose energy very rapidly via 
synchrotron radiation; and photons have a relatively short absorption 
length. 

• What is the acceleration mechanism,? 

The conventional (Fermi) mechanism is the stochastic shock-wave 
acceleration in magnetized clouds. Powerful enough accelerators [3] 
are pulsars, AGN cores and jets or radio-loud quasars. Gamma ray 
bursters may accelerate protons up to lO^'^ eV but they are too dis- 
tant objects [4]. Alternatively, "top-down" scenarios [5] have been 
proposed, based on the decay or annihilation of super-heavy particles 
produced by topological defects [6] or supermassive cosmological relics 
[7]. Another possibility is tlia assumption of exotic physical laws, such 
as breakdowns of Lorentz invariance [8] or general relativity [9] . 

• Where are the sources? 

According to the conventional acceleration mechanism there are very 
few sources at sight within the GZK distance, insufficient to explain 
the number of observed trans-GZK events. Directional and temporal 
decorrelation with the source due to magnetic fields of the order of nG 
has been proposed as a possible explanation [10]. Distant sources like 
BL Lacertae blazers, outside the GZK volume, have received attention 
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recently, since they may emit very high energy photons cascading into 
secondary ones still above the GZK cutoff while approaching us [1 1] . 

To evade most of the problems above, neutrinos are good candidates: 
they are neutral (unbent by magnetic fields, allowing pairing), stable and 
able to propagate unimpeded to earth from very distant sources: even if 
their cross section was of hadronic size, their mean free path would be 
of the order of 100 Gpc among galaxies and 1 Mpc within a galaxy, so 
they can only be stopped by the atmosphere. It is also suggestive that the 
flux of cosmogenic ncTitrinos matches well the flux of observed trans-GZK 
CRs, since the energy-degrading protons photoproduce pions decaying into 
neutrinos. The main problem is that neutrino cross sections are smaller than 
those of hadronic or electromagnetic interactions by six orders of magnitude. 
Two possible types of solutions have been proposed to the small cross section 
problem of neutrinos: 

- "Z-bursts" [12]. Cosmogenic neutrinos are not the primaries but reso- 
nantly create a local flux of nucleons and photons (within the GZK distance) 
with E > Eqzk via scattering with the cosmic neutrino background. The 
resonance is a Z boson. The necessary neutrino mass range fits well with 
the one inferred from oscillation experiments, but large neutrino fluxes are 
required. Similarly, "gravi-bursts" [13] may proceed via resonant Kaluza- 
Klein excitations of the graviton in theories with large compact extra di- 
mensions or in localized gravity models. 

- Enhance the neutrino cross section at high energies. Neutrinos are 
here the primary particles reaching the atmosphere. This will be our ap- 
proach, in the context of a weakly-coupled string theory assuming extra 
space dimensions at the TeV scale. 

2. The effects of extra dimensions 

Extensions of the Standard Model (SM) with extra dimensions offer new 

ways to accommodate the hierarchies observed in particle physics [14]. A 
very attractive possibility would be to bring the scale of unification with 
gravity from Mpianck ~ 10^^ GeV down to the electroweak scale Mew ~ 1 
TeV. This could result if gravity propagates along a (4-|-n)-dimensional flat 
space with n compact submillimeter dimensions (ADD model [15]) or along 
a (4-|-l)-dimensional slice of anti-deSitter space with a warp factor in the 
metric (RS model [16]). These higher dimensional field theories, however, 
must be considered effective low-energy limits only valid below the mass 
scale of a more fundamental theory. And nowadays, only string theory 
[17] provides a consistent framework for the unification of gravity with the 
standard model. 
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In string theory the massless graviton comes as the zero mode of a 
closed string, whereas the gauge bosons and matter fields are the lightest 
modes of open strings. The string vibration modes, with squared masses 

= nMg (integer n > 1), are the so called string Rcgge (SR) excitations. 
The fundamental string scale is Ms = a'^^^"^ where a' is the string tension. 

At energies where the effects of a higher dimensional graviton are un- 
suppressed one expects the presence of its SR excitations giving an effect of 
the same size. This is necessary in order to avoid the pathologies of spin-2 
field theories, such as the quantum field theory of gravity. 




Fig. 1. String graphs for (a) 0{g) vertex of three open strings, (b) 0{k) vertex of 

two open strings and one closed string, (c) 0{cf') tree-level 2-body scattering of 
open strings (no gravity), (d) C(k^) tree- level 2-body scattering with only gravity, 
(e) 0{g^) loop-level 2-body scattering of open strings, topological equivalent to (d) 
and hence = 0{g'^), meaning that (d) is g'^ smaller than (c). 

Now, as emphasized in [18, 19], the exchange amplitude of a closed string 
has an order suppression versus the exchange of an open string, in string 
perturbation theory. This is illustrated in Fig. 1. In consequence, processes 
that receive sizeable contributions from SR and Kaluza-Klein excitations 
of the graviton and also from SR excitations of the gauge bosons will be 
dominated by the second ones, for s < Mg/a} That is, gravity is subleading 
versus gauge interactions and hence the reactions are dominated at tree 
level by open-string diagrams. This is a generic feature in models of higher 
dimensional gravity embedded in a weakly-coupled string theory. 

^ When the ampUtude ylgauge ~ as/Mg < ^grav ~ a^s'^/M^, namely for as > M|, 
gravity dominates and non-perturbative effects become important. 
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3. Phenomenology of TeV strings 

Cullen, Perelstein and Peskin introduced a TeV-scale string model for 
QED in [19]. It contains electrons and photons at low energies and massive 
SR excitations above the string scale. Their results have been generalized 
in order to obtain string amplitudes for neutrino-quark elastic scattering 
[20], mediated in the SM hy a Z boson in the ^-channel. 

The model results from a simple embedding of the SM interactions into 
Type IIB string theory. It is assumed that the 10-d space of the theory has 6 
dimensions compactified on a torus with common periodicity 2nR, and that 
N coincident D3-branes (4-dimensional hypersurfaces where open strings 
may end) are stretched out in the 4 extended dimensions. Such configu- 
ration has N=4 supersymmetry and U{N) gauge symmetry. One assumes 
that the extra symmetry of the massless string modes can be eliminated 
by an appropriate orbifold projection, resulting an acceptable model with 
(at least) the SM fields. The parameters of this theory would be the string 
scale Ms = a'~^^^ « 1 TeV and the dimensionless gauge coupling constant 
g, unified at Mg. Proposals for splitting these couplings can be found in 
[21]. For more general D-brane models see [22] and references therein. 

A tree-level amplitude of open string states on a D-brane is given [19, 23] 
as a sum of ordered amplitudes multiplied by Chan-Paton traces. For the 
processes under study we have 

^(1,2,3,4) = g'^ S{s,t) F^^^^{s,t,u) Tr[tHW + tHW] 

+ 52 s{t,u) F^^^\t,u,s) Tr:[tHHH^ + tHHH^]. (1) 

In this expression, 

_ r(l - a'sni - a't) 
T{l-a's-a't) 

is basically the Vencziano amplitude [24]; (1,2,3,4) label the external par- 
ticles involved in the process: (i^-l^", w^, u^"*) for instance; the Chan- 
Paton factors t"- are representation matrices of U{N); and F^^'^{s,t,u) is a 
factor depending on the vertex operators for the external states and their 
ordering. In our case all the vertex operators correspond to (massless) Weyl 
spinors of helicity (directed inward) + or — , giving 

F-++-{s,t,u) = -d; 

„ 4- 4- / N .u"^ f u u 

F-+-+{s,t,u) = -4— = 4 - + - 
st \s t 

F--++{s,t,u) = -d. (3) 
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Therefore the amphtude for the process i'lul i^lul reads 

n2 



^ S{s, t) Tu43 + - S{s, u) ri234 + ^ u) ^1324 

t u tu 



(4) 



with Tabcd the Chan-Paton traces. To understand the phenomenological 
consequences of this amphtude let us start with the hmit s,t — >^ 0. Since 
r(l) = 1, we have all the Veneziano factors 5(0,0) = 1. The amplitude 
expresses then the exchange of massless vector modes in the t- and the 
■u-channels. The former would correspond to the Z gauge boson, whereas 
the field exchanged in the u-channel is in the (3, 1) and/or the (3, 3) repre- 
sentations of 5'[/(3)c X SU{2)l and has electric charge Q = —2/3. Wc are 
interested, however, in models that reproduce the SM result at low energies, 
with no massless leptoquarks. We obtain this limit if the Chan-Paton factors 
assigned to and vl are such that T1243 — T1324 = —jq and T1234 = T1324, 
where we have used sin^ 9w = 3/8, implied by gauge-coupling unification. 
In terms of ri234 = —a/10 the amplitude (4) becomes 



A 2 2 
5 



- [(1 + a) S(s, t) - aS(t, u)] + -a [S(s, u) - S(t, u)] 
t u 



(5) 



At low s this amplitude is Ao ~ {2/b)g'^ s /t and corresponds to the exchange 
of a Z boson in the t-channel. The Z is then a massless SR mode that 
acquires its mass Mz only through the Higgs mechanism. We shall neglect 
the corrections of order M^/Mg that may affect the massive SR modes. 

As the energy increases the Veneziano factor S{s, t) gives a series of poles 
(at 1 — a's = 0, —1, —2, ...) and zeroes (at 1 — a's — a't = 0, —1, —2, ...). It 
can be expressed as 



'^^^''^"^ a'i + n-1 (a's-n) (n-1)! ■ 



At s = nM| the amplitude describes the exchange of a collection of reso- 
nances with the same mass and different spin (see below). Away from the 
poles the interference of resonances at different mass levels produces the 
usual soft (Regge) behavior of the string in the ultraviolet. Obviously, these 
resonances are not stable and at one loop will get an imaginary part in 
their propagator. When the total width of a resonance (which grows with 
its mass) is similar to the mass difference with the resonance in the next 
level one cannot see resonances and interference effects dominate also at 
s = nM§. 



7 



Let us first analyze the case with a = in Eq. (5). The ampHtude is 
just A{iyLUL ^lul) = (2/5)(7^(s/t) • S{s,t). Near the pole at s = raM|, 

A ^, 2 ^2 nMl jt/Ml) {t/Ml + 1) • .■. • (t/Mg + n - 1) 
•^"~5^ t {n-l)\{s-nMl) ' 

This amplitude corresponds to the s-channel exchange of massive lepto- 
quarks in the (3, 3) representation of 5C/(3)c x SU{2)l with electric charge 
Q = 2/3. At each pole we have contributions of resonances with a common 
mass ^JnMs but different spin, going from zero to the order of the residue 
P„(t) = An • {s — nM|). In this case the maximum spin at the n level is 
J = n — 1. To separate these contributions we first write the residue in 
terms of the scattering angle 6, with t = — (nMj/2)(l — cos^). Then we 
express Pn{0) as a linear combination of the d-functions (rotation matrix 
elements) : 

Pn{e) = \g-'nMl £a;J<oW- (8) 

^ j=o 

The coefficient a'^ gives the contribution to our amplitude of a leptoquark 
of mass nM| and spin J. For example, at the first SR level we find a 
scalar resonance with a? = 1, at s = 2M| there is a single vector resonance 
with a\ = 1, whereas at s = 3M| there are modes of spin J = 2 {a\ = 3/4) 
and J = {al = 1/4). 

The general case with a 7^ is completely analogous, with resonant 
contributions from the terms proportional to S{s, t) and S{s, u). Taking u = 
— (nM|/2)(l+cos 0) and expressing again the residue in terms of d-functions 
we find the same type of resonances but with different a:^ coefficients: = 
1 + 2a, al = 1, ai = 3(1 + 2a)/4 and = 1(1 + 2a)/4. 

4. Application to the iv— nucleon scattering 

From the resonant amplitude u^ul — ^ — i^lUl we can now obtain 
the partial width T:^ = T(X;l i^lul)'- 

" 407r 2J+1 ■ 

Notice that for a given spin J, the variation with n of gives the running 
of the coupling with the energy. We obtain numerically that the coupling 
of heavier resonances decreases like the power law ~ 1/n. 

The partial width T'^ can be used to obtain the cross section (7^(1//, ti/,) = 
a{uLUL — X^) in the narrow-width approximation (NWA): 

47r2 

aiiuLUL) = (2J + 1) 6(3 - nMl) . (10) 
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At each mass level n there is a tower of resonances of integer spin J from to 
n — 1. We find a sum rule for the production rate an{vLUL) = Zlj f^ni^LUL) 
of any of these resonances: 



Eq. (11) is equivalent (for a = 0) to the production rate of a single resonance 
of mass ^JnMs and coupling (2/5)^^ [25]. This is a very interesting result: 
the coupling of heavier SR modes decreases quadratically with the energy, 
but the number of modes (and the highest spin) at each mass level n grows 
also quadratically making j Q^n constant independent of n. 

In the NWA oiyi^uj}) = Y^n j'^i'^LUL — -'^n ~^ anything) is then 
a{i'LUL) = J2n^n{'^LUL)- In this limit the cross section is proportional to 
a collection of delta functions and thus all interference effects are ignored. 
This is a good approximation as far as the total width of a resonance is 
smaller than the mass difference with the next resonance of same spin. Al- 
though the coupling (and any partial width) decreases with the mass, the 
total width of heavier resonances grows due to the larger number of de- 
cay modes that are kinematically allowed. We estimate that contributions 
to from modes beyond n • (g'^/An) « 1 are a continuum. In this 

regime any cross section goes to zero exponentially at fixed angle {s large, 
t/s fixed) and like a power law at small angles {s large, t fixed) [17]. We ne- 
glect these contributions. We keep only resonant contributions from levels 
n < ricut = 50. Our result depends very mildly on the actual value of Ticut- 

To evaluate the total uN cross section one also needs the elastic am- 
plitudes I'Ld'L, ^lur, ^hdR and J^l^l(_r)- ^{uLdi^ i^Ldh) takes the same 
form as the amplitude in Eq. (5) with the changes (2/5, a) (—3/5, a'). 
The massive resonances exchanged in the s-channel are now an admixture 
of an SU{2)l singlet and a triplet. The singlet contribution is required 
in n-even levels, otherwise an SU{2)l gauge transformation would relate 
the parameters a'^ obtained here with the ones deduced from Eq. (5). In 
n-odd mass levels gauge invariance could be obtained with no singlets for 
a' = —{2 + a)/2>. The cross section crn{i^LdL) can be read from Eq. (11) just 
by changing (2/5, a) (3/5, a'). 

The calculation of amplitudes and cross sections for Qr are completely 
analogous. We obtain 




for n odd 



for n even. 



(11) 



' 2 TT ^ 

- ^ h 8{s - nM^) for n 
for n 



(12) 
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The cross sections (Jnii^LdR), crn{i^LUR) and (Tni'^^LdR) coincide with the 
expression in Eq. (12) with the changes (2/5,6) (1/5,6'), (2/5,6) — 
(2/5, a) and (2/5,6) (3/5, a'), respectively. For the left-handed anti- 
quarks, (Jnii-'LML) and cJn('^Lt^L) Can be read from Eq. (11) by changing 
(2/5, a) ^ (2/5,6) and (2/5, a) ^ (1/5,6'), respectively. 

Now the total neutrino-nucleon cross section due to the exchange of SR 
excitations can be very easily evaluated. In terms of parton distribution 
functions q{x,Q) {q = qL,R,QL,R) ^ nucleon (N = {n+p)/2) and the 
fraction of longitudinal momentum x, it is 

<-ln) = y.i:-^-<im), (13) 

n=l Q S 

where x = nM^/s, s = 2MNEy, = nMg and dniyhl) is the factor 
multiplying the delta function in the cross section (Jn{i^LQ)- 

10-29 

10-31 
„ 10-32 

a 

- 10-33 
b 

10-34 

10-35 
10-36 

10-37 

10^ 103 lO'' 10^ 10" 10^ 10*^ 10^ 101° 10" 1012 10" 

[GeV] 

Fig. 2. Neutrino-nucleon cross section versus the incident neutrino energy E^. The 
SM contribution (solid) includes neutral and charged current interactions. The SR 
contribution is for Ms = 0.5 TcV (dashes) and Ms = 2 TeV (dots) for the cases 
(i) a = a' = 6 = 6' = and (ii) a = a' = 6 = 6' = 5. 

In Fig. 2 we plot the neutrino-nucleon cross section at energies from 
102 to 1013 GeV for Ms = 0.5, 2 TeV. We have used the CTEQ5 parton 
distributions in the DIS scheme [26] extended to x < 10~^ with the methods 
in [27]. We show the SM cross section and the string corrections for a = 
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a! = h = h' equal and 5 (notice that in tlie first case there are no s-channel 
resonances mediating the vlQr amplitude). The modes beyond ricut — 

50 

are not included, since there we expect that the narrow width approximation 
is poor. 

As a final comment, let us remark that the SM neutrino- nucleon interac- 
tion probes values of x much below the current HERA DIS data (a:;^^!^^ 
10~^) for ultrahigh-energy neutrinos, 



(TSM 167r -y/s \GeV 

This is not the case for the process mediated by SR resonances, where the 
probed values of x are above 

taking Ms = 1 TeV. 



5. Conclusions 

Cosmic rays hit the nucleons in the atmosphere with energies of up to 
10^^ GeV. If the string scale is in the TeV range, these cosmic rays have the 
energy required to explore the fundamental theory and its interactions. In 
particular, ultrahigh-energy neutrinos are interesting since they can travel 
long distances without losing a significant fraction of energy and arc not 
deflected by magnetic fields. In addition, the SM interactions of a neutrino 
are much weaker than those of a quark or a charged lepton, which makes 
easier to see deviations due to new physics. 

With this motivation we have analyzed the string uN cross section at 
energies much larger than the fundamental scale Ms ~ 1 TeV. In a weakly- 
coupled string theory the process is given at tree level by open-string graphs, 
whereas gravity effects appear as a one-loop correction.^ We have fixed the 
arbitrary parameters of the model imposing phenomenological constraints, 
namely, the massless SR modes must account for the standard model par- 
ticles. Four Chan-Paton traces remain as free parameters. The massive 
SR modes include leptoquarks that resonantly mediate the process in the 
s-channel for energies above the string scale. The presence of massive lep- 
toquarks is not a peculiarity of our toy model but a generic feature of any 
string model, and has to do with (s,t) and/or {t,u) duality of the open- 
string amplitudes (see e.g. Eq. (4)). 

^ Gravity effects in the context of weak-scale string theories at high energies have been 
considered in the Regge picture by [28]. 
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A very simple sum rule for the production rate of all the s-channel 
Icptoquarks, with different spins in the same mass level n, makes possible the 
calculation of the total uN cross section in the narrow-width approximation. 
The effect of these leptoquarks is not just a correction of order M^/Mj to 
the SM cross section, as one would expect on dimensional grounds. Such SR 
excitations give a contribution that can dominate for Ms ~ 1 TeV. However, 
for the expected flux of ultrahigh-cncrgy neutrinos it seems unlikely that the 
cosmic ray events observed above the GZK limit correspond to the decay of 
string resonances produced in vN scattering. A similar conclusion has been 
recently drawn in [29] for graviton-mediated vN scattering and black hole 
production in TcV-gravity models. 

Nevertheless, since neutrinos are very penetrating particles, the enhance- 
ment of the cross-section may make possible the detection of horizontal air 
showers produced by ultrahigh-cncrgy CRs in upcoming experiments [30]. 
Furthermore, in second-generation neutrino telescopes, able to detect neu- 
trinos in the TeV to PeV range, there is also a chance to probe this and 
other models of TeV-scale quantum gravity at more moderate energies [31]. 

It is my pleasure to thank the organizers of the XXV International School 
of Theoretical Physics for their kind hospitality and the stimulating and 
enjoyable atmosphere of the event. The author acknowledges F. Cornet and 
M. Masip for carefully reading the manuscript and a fruitful collaboration. 
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